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Abstract 

The structure of counterterms in higher derivative quantum gravity is reex- 
amined. Nontrivial dependence of charges on the gauge and parametrization is 
estabhshed. Exphcit calculations of two-loop contributions are carried out with 
the help of the generalized renormgroup method demonstrating consistency of 
the results obtained. 

1 .Introduction 

As well known, not all of the problems of the quantum field theory are exhausted 
by the construction of S-matrix. Investigation of evolution of the Universe, behavior 
of quarks in quantum chromodynamics etc. require the introduction of more general 
object — the so called effective action. Besides that, the program of renormalization 
of the S-matrix itself has not yet been carried out in terms of the S-matrix alone. 
Renormalization of the Green functions, therefore, is the central point of the whole 
theory. Given these functions one can obtain the S-matrix elements with the help of the 
reduction formulas. In this respect those properties of the generating functionals which 
remain valid after the transition to the S-matrix is made are of special importance. 

We mean first of all the properties of the so called "essential" coupling constants in 
the sense of S.Weinberg They are defined as those independent from any redefinition 
of the fields. In the context of the quantum theory one can say that the renormalization 
of "essential" charges is independent from renormalizations of the fields. Separation of 
quantities into "essential" and "inessential" ones is convenient and we use it below. 

In this paper we shall consider the problem of gauge and parametrization depen- 
dence of the effective action of i?^-gravity. 

There are two general and powerful methods of investigation of gauge dependence 
in quantum field theory. The first of them uses the Batalin-Vilkovisky formalism 
P, 1^, ^ and is based on the fact that any change of gauge condition can be presented 
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as a (local) canonical (in the sense of " antibrackets" [Q) transformation of the effective 
action. This canonical transformation induces corresponding renormalized canonical 
transformation of the renormalized effective action. This leads to the following result: 
the renormalization boils down to the redefinition of the coupling constants (which 
are the coefficients of independent gauge invariant structures entering the Lagrangian) 
and some canonical transformation of the fields and sources of BRST-transformations. 
The second approach ^, || consists in the introduction of some additional anti- 
commuting source to the effective action in such a way that the Slavnov identities 
for the corresponding generating functional of proper vertices connect its derivatives 
with respect to the gauge fixing parameter and to the mean fields (and sources of 
BRST-transformations) . 

The second method was used in to prove the gauge independence of the gauge 
invariant divergent parts of the effective action up to the terms proportional to the 
classical equations of motion of the gravitational field. Together with the general result 
of the first approach mentioned above this would imply some far-reaching consequences 
concerning renormalization of the fields. For example, one could conclude that the 
canonical transformation corresponding to a change of the gauge condition should not 
be renormalized. Unfortunately, this is not the case. We will show in this paper that 
the aforesaid result of |Q holds at the one-loop level only. Introduction of the additional 
source mentioned above requires also introducing of some additional terms needed to 
compensate divergences which arise because of the presence of the new source. As a 
result corresponding Slavnov identities impose only some constraints on the divergent 
structures from which a nontrivial dependence^ on the gauge follows already at the 
two- loop level. 

Our paper is organized as follows. In sec. 2 we determine possible divergent struc- 
tures which are originated due to the presence of the new source and obtain correct 
Slavnov identities. In sec. 3 we calculate explicitly the divergent part of the effec- 
tive action at the one loop level in arbitrary (linear) gauge and the special class of 
parametrizations. In sec. 4 we calculate the divergent as ^ {e being the dimensional 
regulatoij^) part at the two-loop level with the help of the generalized renormgroup 
method and show that the results obtained in sees. 3,4 satisfy the relations derived in 
sec. 2. 

We use highly condensed notations of DeWitt throughout this paper. Also left 
derivatives with respect to anticommuting variables are used. The dimensional regu- 
larization of all divergent quantities is supposed. 

2. Generating functionals and Slavnov identities 

a. Action, gauge fixing and parametrization 

Let us consider higher derivative quantum gravity described by an action which 
includes the minimal set of terms added to the usual ones of Einstein to ensure the 
power counting renormalizability of the theory^ 

^I.e. dependence which can not be presented as proportional to the equations of motion. 

^We set 2e = d — 4, d being the dimensionahty of the space-time. 

^Our notations are R^i, = = da^'jlu - R = RuiyQ^" , 9tiu = sign{+, -). 
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5-0 = 1 d^x V^iaiR^ + (5R^,R^' ~ ~ ^A)), (1) 

where ai and /3 are arbitrary constants satisfying only /3 7^ 0, 3ai + /3 7^ which 
imply that the graviton propagator behaves like for large momenta (see [^); k is 
the gravitational constant and A is the cosmological term. 

The corresponding equations of motion are 



-2/5i?^^i?^"'^^ - (2ai + URg^l^ - pUR^I^ 

+ (2ai + m^^" - ^Rg'"' + ^Ag'^" + ^R'^" = 0. (2) 

Renormalizability of this theory was proved in in the case of the so called un- 
weighted (or weighted with a functional containing fourth or higher derivatives) har- 
monic gauge condition. The proof in more general case boils down to the proof of the 
so called locality hypotheses. In [|TU[ its validity was shown most generally. 



For our purposes it is sufficient to consider the harmonic gaugeQ following Stelle 



= FfKp = d'h^, = 0, (3) 

where hfj^^, denotes some set of dynamical variables describing the gravitational field. 
We recall that in the theory of gravity a natural ambiguity in the choice of such a set 
exists because the generators D"^, of gauge transformations of variables constructed 
from the metric g^u (or g^"^) and its determinant g = detg^v in any combination have 
simple form linear in fields and their derivatives. For general constructions of this 
section it doesn't matter what choice we make. We only note that the gauge in the 
form of d^) will always correspond to the set of dynamical variables]^ which enter the 
so called reduced expression of the metric expansion (see sec. 3a, eq.(p^). Thus the 
BRST-transformations of the Faddeev-Popov effective action with the gauge fixing 
term being — ^F^DFq,, expressed in terms of these standard variables are 



5C^ = -A-^DF^X, (4) 
where A is an anticommuting constant parameter. 

b. Green functions 



We write the generating functional of Green functions in the extended form of 
Zinn- Justin modified by Kluberg-Stern and Zuber 

^We use the flat-space metric tensor 77'^" = diag{+l, —1, —1, —1) to raise Lorentz indices. 
^They will be referred to as standard variables. 
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J dh^,dC^dC^exp{i{t{h^,, a, C\ K''^ L'^) 

+YF„C'' + P^'C^ + C^Pr + T^^h^u)} (5) 

where 

K^"{x) (anticommuting), L"{x) (commuting) are the BRST-transformation sources 
and Y is a constant anticommuting parameter. 

Let us first consider the structure of divergencies which correspond to the extra 
source Y. Power counting gives for the degree of divergence D of an arbitrary diagram 

D = 4:-2n2-2nK -nL-2nY - Ec -2Ec, (6) 

where n2 = number of graviton vertices with two derivatives, nK,L,Y = numbers 
of K,L,Y vertices respectively, Ec and Eq = numbers of external ghost and antighost 
lines. 

Also from the expression (|^) we see that one can ascribe the following ghost numbers 
Ng to all the fields and sources: 

Ng[h]=0,Ng[C] = +l,Ng[C] = -l, 

Ng[K] = -1, Ng[L] = -2, Ng[Y] = +1. (7) 

Now from and one can see that there are three types of divergent structures 
involving Y- vertex: Y K, Y C* and Y L C, each of which may have arbitrary number 
of external graviton lines. As far as we have adopted the standard covariant approach 
thus only Lorentz covariant quantities may appear and therefore we have for the general 
form of the above structures 

YK'^'^P^,, (8) 



YC-'d^Q,, (9) 

and 

YL^'C^M^r, (10) 

where P, Q and M are some Lorentz-covariant tensors depending on h^i, alone. 
Thus to renormalize the Green functions we must introduce corresponding coun- 
terterms and consider the new generating functional 
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+YL''C^M,r + ra + C^Pr + T'^'^ V)} (11) 



instead of 



c. Slavnov identities 

Let us proceed to successive renormalization of Green functions corresponding to 
(|TI|). We will first consider the case when (7*^ = g^^ are chosen as a parametrization 
of the gravitational field. Then the general result will be clear. 

To ensure renormalizability we work with a generating functional (pH]) from the 
very beginning. We will see below that Slavnov identities determine the structure of 
the polynomials P and Q completely. They turn out to be 

PfMu = a{ri^^ + hf,^), (12) 

Qf,u = ah^^, (13) 
a being some divergent constant. Thus we set 



at the zero order. Then inclusion of the counterterms (12,13) is just a multiplicative 
renormalization of the source Y. 

a. One-loop order 

To obtain Slavnov identities at this order we perform a BRST-shift (^) of integration 
variables in (pH]) 



5MW _..„„\ / 6 



+^rC"^Fri^^,a]ea;p{2(S + YF^C'' + YK^'^ir]^, + h^,) + FL'^C^M^^) 
+P"C, + C^l3r + T^'h^,)} = 0. (14) 

Our aim is to find the A-dependence of the gauge invariant terms only. Terms 
containing M'^) in (|l^) depending on anticommuting fields C„ and source L'^ are unim- 



portant in this respect and we replace them simply by "+..." in what follows because 
these terms will be omitted in the end of the calculation anyway. 
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Using the ghost equation of motion 



exp{t{t + YF^C' + YK'^'^ir]^, + V) + - 
introducing the generating functional of proper vertices F 



f[h^,,C^,C^,K^\L'',Y] 



w 



-i InZ, 



and noting that 



h 



_ sw 

6W 



dV dW 



dA dA ' 

we rewrite (p^ as the Slavnov identity for F 



6r 



Sh 



1 



6T 



+ Y{r]^, + h 



+ 



A "5a dA 5a 



0. 



To simphfy eq. ([l9|) we introduce the reduced generating functional 



F = f + ^i^aDF" - YK^\^^, + V) - YF,C\ 



Then eq. (|l^) reduces to 



5F 5F 5F 5F 

+ 



bh^^ bK^'^ bC„ bL'' 

dA bKi"' 



0. 



The ghost equation of motion written in terms of F is 
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- ^ = 0. (21) 

Now let us separate Y-independent part of F from the part linear in Y 



r = ri + yr2. (22) 

Then eq.(pO|) gives an ordinary Slavnov identity for Fi 

i^ii^^ + ^^ = (23) 
and an equation involving r2 



bh^^bK^'^ bh^^bK^'^ bC^bL" 



■2A^ + A-i^ + .„^0, (24) 



bC^bL'' dA bK^''' 

Finally, we omit all but the terms depending on h/j^^, only and obtain in the first 
order 



' +2A— -— + ... = 0, (25) 



bhf,^ bKt"" dA 

where fl denotes the gauge invariant part of Fi and the superscript div{l) denotes 
the one-loop diverg;ent part of the corresponding quantities. 

As we know Tg™^^^ = i^'^^P^i^ P^^^ being some divergent polynom in /i^,^. 

Thus, dropping out the terms proportional to K'^'^ again and the symbol " + ... " 
we obtain the following equation for the gauge invariant terms 

Qdiv{i) ^Yie effective 

action^ 

The left hand side of this equation is gauge invariant thus so is the right hand side. 
Therefore P^^^ has the form mentioned above. 

To make the Green functions finite at the one-loop level we must redefine the initial 
effective action E 



S ^ S(i) = S - rf (27) 

^We will see in sees. 3, 4 that the non-gauge-invariant terms in Tf™ depending on /i^^ only are 
absent. 
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and the source YQ 



r^r(l-a«). (28) 

As explained in ^ subtraction of rf^^^^ boils down to a redefinition of all the fields 
in such a way that T.^^^ is invariant under renormalized set of BRST-transformations 
for which we do not introduce new notation. 

p. Two-loop order 

We perform a renormalized BRST-transformation of integration variables in the 
generating functional of Green functions finite at one-loop level 



^[1] ^j^u^ pa^ j^^u^ j^a^ Y] = 

J dh^^dC„dC^exp{i{t^^^ + F(l - a^^^)F„C'' + (1 - a^^^)Y K^'^ [r]^^ + V) 
+... + ^'^a + C^f3r + T^'^'Ku)} (29) 
and obtain the following Slavnov identity 

j dh^,dC,dC^[(T^'' + Y{1 - a^'y)K>^'') + Y{1 - a^'^){r]^, + V)) 



^ + 4 ;5^nF-,M-^ _ 2Y(1 - a(^))A f + ^ 



div{l) ' 



^ ' A'"' ' 5T^''' ' ' UA ' rfA 

+iY{l - a^^^)C''F^''D';:,Ca + ...]exp{2(S« + Y{1 - a^^^)F,C'' 
+r (1 - a(i))ir^^(r/^, + V) + - + + C*^/?. + T'^'' V)} = 0. (30) 
To evaluate the termQ 

j dh^^dC^dC^YA exp{i(S(^^ + r(l - a(^))F,C" 

+ (1 - a(^)) W^(r/^, + V) + ... + ra + ^"/3r + T^'^ V)} (31) 
we use eq.(^) and equation of motion of h- field which is obtained from (^9l)P| 



Y J dh^^dC^dC 



5rf 



exj9{i(S(^) + F(l - a(^))F,C"^ 

+y(l - a(i))ir^'^(r/^, + V) + ... + P^C^ + C^Pr + T^^'K.)} = 0. (32) 



^We should also include counterterms of the type YLC, but they are irrelevant to the issue and 
replaced by " + ..." as we have mentioned above. 

Again evaluation of the gauge invariant part of A — ^ — is needed only. 
^We use the property = 0. 
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In the two-loop approximation we may write 



^ exp{i...} = ct^^^^T" 

Finally, using the ghost equation of motion 



/ dh.,dC^dC^--^ exp{i...} = a(^)^-i ZW. (33) 

dh„,, dh,,u 



F^T^ - ^ = 0, (34) 

written in terms of the one-loop finite reduced generating functional of proper 
vertices 

rW = fW + ^^^n^" - YK^\r]^u + Ku) - YF,C\ 

we rewrite the rest of the eq.(|5DD as in section a and obtain the following Slavnov 
identity for one-loop finite proper vertices, valid up to two-loop order 

where terms explicitly dependent on K^^ and (7^ are included in " -|- ... " for 
simplicity. 

Again the separation 

rw = + yr™ (36) 

gives an ordinary Slavnov identity 

and an identity involving Fg^' 

5Tf 5fW ^ 5rf 5rf ^ 5fW 5fW ^F™ ^F™ 



+2A— ^ + a«^-^(r7^, + V) + - = 0- (38) 

Thus for two-loop gauge invariant divergent part fi'^l'^™*^^) of the one-loop finite 
generating functional of proper vertices we have 



5h^, ^'^ dA 
+a(') —7 (r/^. + V) = 0. (39) 
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Again it follows from eq.(g|) that Pj^) = a^^\r]^^ + h^^). 

Thus we obtain the following identity for one- and two-loop divergent gauge invari- 
ant parts of the effective action 

2A = -«^^^^(^/^- + V) + « ^(^M^ + V)- (40) 

Had we used any other parametrization of the gravitational field, P'-^-* and P*^^-* 
would have such a form that provides the gauge invariance of the product j^P^^, 
where h^^j denotes the set of standard variables. Therefore result (|40D holds in general 
if h^y denotes a quantum part of the covariant components of the metric field. 

Thus we see that in the presence of the new source Y the renormalization procedure 
differs from the usual one substantially. Although the renormalized Green functions 
satisfy the same Slavnov identities as the bare ones, the renormalization equation (of 
the type (|40|) ) for their divergent parts in (n+l)-th loop order can not be obtained 
by simple omitting of the finite parts of the Slavnov identities for the Green functions 
renormalized up to n-th loop order. The correct procedure presented above leads to 
the Slavnov identities which just impose some nontrivial constraints on the form of 
gauge-dependent divergent structures of the Green functions. 



3. Calculation of the one-loop divergent part of 



In previous section we have obtained the relation ( PD| ) which identifies [modulo terms 
proportional to the equations of motion of h field) the A-derivatives of the two-loop 
gauge-invariant divergent part of the effective action with the variational derivatives of 
the corresponding one-loop part up to some coefficient being defined by divergent parts 
of diagrams with one insertion of the Y vertex. To prove this coefficient is not zero we 
present explicit calculation of the values r^™''^'* and pW"^*^*^^) arbitrary gauge of the 
type (§) and arbitrary parametrization with the only restriction of linearity of group 
generators. We prefer this way to direct computation of diagrams with Y insertion 
because it allows to verify the relation (^). 

a. Arbitrary parametrizations 

In general the metric is an arbitrary function of dynamical variables. The only 
restriction is that this function must be nondegenerate. For example, if dynamical 
variables are chosen as g*^^^ = g^iu^—gY.g = detg^^ then we should avoid the case of 
r = — i; otherwise detg*^ = 1 and one more independent variable must be introduced 
in addition to the set of g 



To calculate one-loop divergences the background field method is used |]TT], 0, 
Accordingly, we should first find expansions of all the quantities entering (0) in powers 
of dynamical quantum variables h^i, around the background field g^^u up to the second 
order. 

Now we note that the form of the graviton propagator is, of course, parametrization- 
dependent and this dependence complicates all calculations considerably. However, it 
is fictitious in the sense that always can be removed by linear redefinition of quantum 
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variables. Such a change doesn't mix different orders in the loop expansion and there- 
fore doesn't alter the values of the one-loop divergent part in particulaiQ. We will 
show below that our calculations are highly simplified if the linear part of the metric 
expansion is chosen to have the simplest form /i^j, 

g^^^ = Qfiv + V + ahh^^ + hh^.^h'^ + ch^g^^ + dh^ph'^'^g^^ + 0{h^), (41) 

where a, b, c, d are arbitrary constants; g^^ denotes the full metric field; all raising of 
indices is done by means of the inverse background metric g'^": g^ag°"^ = ^'jl'i h = h^^g^^. 
Any parametrization g* somehow constructed from the metric g^y and its determinant 



has a background expansion reducible to pTf). 



To show an advantage of such a choice of the background metric expansion we 
note that it is only the linear part of this expansion which in fact contributes to the 
curvature R'^ap expansion. Really, this tensor has the following structure 

R = dV- dV + TT-TT. (42) 

Suppose at the moment that we have chosen our coordinate system in such a way 
that r = at any fixed point of space-time.|^ Then 

R = dT-dT + riFi - riFi + d{Ti + T2) - 9(ri + T2) + oih"), (43) 

where subscripts 1 and 2 denote parts of F of the corresponding powers in h. We 
may rewrite this expression in explicitly covariant form as 

R = R + V(ri + Ta) - V(ri + Ta) + T^T^ - T^T^ + ©(/i^), 12 (44) 

valid therefore in every coordinate system. In terms of the Lagrangian linear in 
curvature scalar all the full derivatives of the second order may be dropped out. In 
quadratic terms these derivatives are multiplied by the zeroth order quantities R^v, R 
etc., when the second variation of the action is being calculated. Integrating by parts 
one can easily verify by means of power counting that these terms do not contribute 
to the one-loop divergent part of the effective action. 

It follows from the above discussion that the dependence on the parametrization 
appears only in terms without h-derivatives in the second variation of the action if the 
reduced expansion (^) is used. 

b. One-loop invariance on-shell. 

To calculate one-loop divergent part of the effective action in arbitrary gauge and 
parametrization we shall use the fact that the dependence on parameters A, a, 6, c, d 
appears in the terms proportional to the equations of motion only. As far as the 
A-dependence is concerned the corresponding result follows directly from eq.(^). 

To prove the on-shell independence from a, 6, c, d we note first of all that these 
parameters appear in the second variation of the action only in terms having the form 



^''The corresponding Jacobian is exp{6{Q).. ) = 1 in the dimensional regularization. 

Recall that F and T are constructed from g and g^i, respectively. 
-"^^Note that Fi ^ are tensors. 
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where [g^u]2 denotes the second order part of the reduced metric expansion (^Tj). 

Next, calculating generators of the gauge transformations of dynamical variables 
belonging to the class of parametrizations described above and passing to the set of 
standard variables again one easily sees that these generators just coincide with the 
ordinary ones of the metric field transformations, i.e. they are a, b,c,d - independent 
and therefore so is the ghost contribution. 

Thus the on-shell invariance is proved. 

c. Background field method. 

According to the background field method we separate the quantum field part h*^^ 
liTom the external field 5'*^, 

Then we expand the metric field g^i, in powers of /i*^ and perform a linear trans- 
formation on /i*^ bringing this expansion to the form of (|4l|). 



Imposing the background Lorentz gauge on the quantum field h 



111/ 



FM = Ff{g)Kp ^ V' (45) 
we have for the generating functional of Green functions 



1 



exp{t{S,{g, h) - F"{g)V'F^{g)V^ + C^Fr{g)Dl^C^ + T'^^V)}- (46) 



2A 

We suppose that the background field g^j^y — rj^j^u and the source T^^^ are absent out 
of some finite region of space-time. Integration is carried out in all fields h^^ tending to 
zero at infinity. We do not introduce background ghost fields nor their sources because 
renormalization of these fields plays no role in this section nor in sec. 4. 

In the one-loop approximation we expand the gauge fixed action 

S,s ^ So{g, h) - ^F-[g)V^F^{g)^ 
around the extremal h satisfying the classical equations of motion 

5Sgf{g,h) 



T^^ = (47) 



up to the second order and obtain 
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Z[T^^] = exp{i{Sgf{g, h) + T^^^^h^} [det g,,V^f 

J dh^,dC„dC-detFl''{g)D-^^{g, /i)ea;p - Ku){Kp - Kp)^^^) 

As far as we have supposed the background field g — r] and the source T to disappear 
out of some finite region of space-time one can choose a solution h of eq. ( ^Tf ) to be zero 
at infinity. Thus the shift of integration variables h h + h doesn't change boundary 
conditions for h and we have for the generating functional of connected Green functions 



W^-t InZ = S,f{g, h) + T^^ V + ^ ' 

I Tr In Fr{g)D'^,,{g, h)-'-Tr In g^,V\ (49) 



To perform a Legendre transformation we calculate 



_ _ 

and 



T{g, h) = W{g, h) - h.^T^'' = Sgf{g, h) + Tr iJ'^'^^^' 



2 5h,,,,Sh. 



« ~ - -2 



t Tr In F!^'{g)D;,{g, h) - -Tr In g,,W 

^ Tr /n^^^^T^ - ^ Tr /n Tr(^7)I^^.(^?, /^) - ^ Tr /n ^7,.V^; (51) 
the eg. (^71) was used in the last passage. 

Obtaining the relation (^) in sec. 2 we used the fiat background r^^;^. Had we started 
with an arbitrary background metric g^i_y instead of ?7^jy we would have modulo terms 
proportional to the equations of motion 



^^^^ = 



We wrote g^.^ instead of g in ( ^2]) because it is sufficient to verify this relation in 
case h^y = 0. 

d. Calculation of 
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Let us first reveal some "essential" properties of charges. 

As i?^-gravity is renormalizable we can write the fi'^™^^^ in the form 



1 r 

^ j d^x (ciR^^R^^, + C2R^ + c-iR + C4A + cg 



(53) 



Qdiv{l) _ 

327r% 

where Cj, z = 1, 5 are some gauge and parametrization dependent coefficients. 

As we know from sec.b. fi'^™*^^) is gauge and parametrization independent on shell. 
It is obvious that the only scalar which can be constructed from equations (||) to 
transform Vt'^'^'^^ [W^ is 



1 



{R - 4A) = -2(3ai + (3)V'^R. 



(54) 



It follows from these simple facts that Q,i = 1,2,5 and the combination 4c3 + C4 
do not depend on A, a, fe, c, d. 

Thus we may simplify the calculation of 1]*^*^^^ in arbitrary gauge and parametriza- 
tion if divide it into two parts: 

1. Calculation of n'^*''^^) in the case of the simplest gauge and parametrization. We 
choose g*^^ = g^^, and the minimal gauge. 

2. Calculation of the coefficient C4 alone in arbitrary gauge and parametrization. 
In this part we may obviously consider the space-time as fiat. 



The correct result of the first part of our program was obtained in |T4 
1 



where 



327r% 



J d^x (ciR^'^R^^ + C2R^ + o^R + C4A + cg) , 



(55) 



133 lOaf lOai 291 

Cl = -TTTi C2 = + ^77, —1 C3 



1 



10 



1 



C4 



(3ai + /5)P 



6/5 60 
28ai 



/5 



+ 9 



,C5 



(3ai + /3)P 
3 



30af 53ai 21 



+ 



2/? ^ 4 



/?2 

15a? 5ai 7 
H + - 

2(3^ P 8 



(3ai + /3)2fc4 

Calculation of C4 in the flat space-time is presented in Appendix A. Combination 
of the two results gives 



1 



where 



327r% 



d'^x (ciR'^'R^^ + C2R^ + C3R + C4A + C5 , 



(56) 



133 



10 2 5 773 



10 ' 9 " 3" 180' (3^k^ V2 ' 8a2 



C4 



4C3 + C4 



3q; 









K 


'2A 


a/ 


40a 




3a 




, a 




3 





i; ( 14A + - + 20 

a 



+ l,M = a-F4c+^,f = 6-F4(i--. 

4 2 
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4. Calculation of the two-loop divergent part of 



As follows from eq.(PDD the nontrivial dependence on the gauge parameter A (i.e. 
which is not zero modulo equations of motion) is contained in terms proportional to 
^. To calculate the latter we use the renormgroup method. It is very convenient to 



apply the generalized version of the renormgroup equations given in WB, |T6|. For the 



sake of completeness we give an account of this method following [|l| 



a. Generalized renormgroup method 

The idea of this approach is to obtain renormgroup equations without explicit 
distinguishing of different charges, i.e. in terms of the whole Lagrangian. 

Let us consider the bare Lagrangian as a functional of the initial Lagrangian L 



(57) 



n=l 



where symbol An{L) means that the corresponding counterpart is calculated for 
the Lagrangian L. Independence of the /^from the mass scale implies 



/3(L)= L— -lUi(L), 



(58) 



(^L^-lj A„(L)=/3(L)^A„_i(L), 

where the so called generalized /3-function /3(L) is defined by 

o dL 



~eL + (3{L), 



Lb 



(59) 



We don't have to muse upon the concrete sense which the operation ^ possesses. 
Using the loop expansion of An 

oo 

A^{L) = ^ A„fc(L), 

k=n 

and noting the homogeneity of functionals Ank{L) 

AnkiXL) = Ai-'A„fc(L), 

(A being a constant) 

we can express the operations Lj^ and P{L)j^ in terms of the ordinary differenti- 
ation 



6 d 
L—Ank{L) = —AnkiXL) 
oL oX 



'l-k)Ank{L) 



(60) 



A=l 
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and 



r 7 

/3(L)— A„fc(L) = —Ank{L + x(3{L)) 

whatever meaning has to be assigned to 
Thus we obtain 



x=0 



(61) 



A=l k=l 



(62) 



'd_ 
dX 



- 1 A„,(AL) 



k=n 



A=l 



k=n—l 



1=1 



To relate 74„„ and we substitute L — > .^^^L in 

respect to ^ n — 1 times and set = 0. 
The result is 



(63) 

x=0 

differentiate with 



nAnn{L) = -^An-l,n~liL + xAu{L)) 



(64) 



x=0 



b. Calculation of n^^jf^^^ 



To apply the relation 



to the case of 



A22{L) 



1 d 

2 dx 



Aii{L + xAii{L)) 



(65) 



Leff = v^(«ii?' + f^R^uR'"'' -^(R- 2A)) 



(66) 



we note first of all that the gauge-fixing term is not renormalized if the linear gauge 
is used (see e.g. [||, |l3])- Also the renormalization of the ghost part of the effective 
action is immaterial as long as only the one-loop expression is needed in 

Thus to calculate Q^y^r*'^'' we rewrite Lgf + xAu (L) as 



^|(ai + xc2)R' + {[3 + xci)R^,R'''' 



R-2 



, 1 , _i / A XC4A + xc^ 
^F-"^^)" If + 2 



(67) 
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apply the one-loop result (|^), differentiate it with respect to x and set x = 0. 
The result of this calculation is presented in Appendix B. 
Now we are in position to verify the identity (|52[). 

As follows from the result (173|) on the mass-shell the left hand side of fp2D is 



A = {2A^w^ + Aw{-20va-^ + 20v + Swa'^ + 40/3a - l/3a"^ - 10)} , 



Aw , , o 
^ = -^(V4a'' + 5), = n + Iv, 

while (^Bp gives for the right hand side 



327r% 



J (fx (iA + B) , 



(69) 



i = -L |2Aw; - 20fa"^ + 2Qv + 3wa"^ + 40/3a - l/3«"^ - loj 



B = -g4;{l/to-H-5} 



We see that the eq.(^) is really satisfied and the coefficient a*^^^ turns out to be 
equal to ^^r^- Note that '^"'^j^"'^^^ is not zero even if the unweighted gauge condition 
A ^ is used. 



Conclusion 



We have shown in this paper that generally the divergent parts of the effective action 
of i?^-gravity depend on the gauge and parametrization nontrivially — this dependence 
can not be presented as proportional to the equations of motion. The renormaliza- 
tion procedure in the presence of the new anticommuting source Y turned out to be 
more complicated then the usual one: the renormalization equation corresponding to 
the modified generating functional can not be obtained by naive extracting of diver- 
gent terms in Slavnov identities. We have considered the renormalization of modified 
Green functions at one- and two-loop levels and obtained renormalization equations 
corresponding to the insertion of Y-source (eq. (p6| , ^0D ). Also explicit calculation of 
the one- and two-loop divergent parts has been carried out confirming our results and 
demonstrating that the nontrivial gauge dependence of the divergent parts of the effec- 
tive action actually exists in arbitrary (Lorentz) gauge and arbitrary parametrizations 
except those satisfying w = CQ. 

^•^Construction of the parametrization satisfying w = see in [l9| . 
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We emphasize that this nontrivial dependence is due to the presence of the Einstein 
term in the Lagrangian. Had we considered a theory with the Lagrangian containing 
the higher derivative terms only we wouldn't have had such a dependence. 

Our conclusion does not contradict the equivalence theorem pO[ in view of the 



general results of [0, pT|. Their validity in the present case is verified in appendix B. 



However, these results do not allow to say that the renormalization of the coupling 
constants is independent from the renormalization of fields (as in the case of two- 
dimensional chiral theories for example), because renormalization of the Newtonian 



gravitational constant k can not be separated from renormalization of the gravitational 
field: one can always perform additional redefinitions of the constant k and the metric 
field which compensate each other. This is a consequence of the fact that k is an 
"inessential" coupling constant. 

Finally, we note that our results are in agreement with the general statements of 
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Appendix A 

In this appendix we present calculation of the one-loop divergent part of the effective 
action in the fiat space-time. 

According to algorithm derived in Jl^ we should first calculate the part of S'^Sgf 
with four derivatives 



9f 



= ha^h (ai + - + ^1 + -h.^ah"' + A^aA" ( - - ^ 
4 I ^ 4 2A / 4 I 2 2A , 



Then we substitute — n^, being a vector with = 1, and calculate the 
"propagator" (Kn^^)"^'^^ which is the inverse of the operator {Kn)^y^af3 = sh^^sh p • 

( j^^-i)a/3,75 = I /2{g'^^gl^^ + g'^^g^^) + Aig'^'^g'^^ + Bi{g''^n'^n^ + g'^^n^n^ + g'^^n'^n^ + 
gl^^n^rC) + Ci{g"^^n'n^ + g'^^rfn^) + Dirfn^rCn^ , where 

Ai = -{A + AAB + AD- C^) /Z, 

5i = -5/(l + 25), 
Ci = {AAB + AD-C- C^)/Z, 
Di = -{1QAB + AAD + D + AB-AC^)/Z + AB/{1 + 2B), 
Z = l + AA + AB + 2C + D + ?,{AAB + AD - C^), 
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coefficients A, B,C, D being defined from L 



gf- 



4 , 4ai 2 



C 



1 



D = 2 + 



P /3A' 
4ai 



We have multiplied the initial Lagrangian by ^ for convenience. 

Second, we calculate the part W of containing two derivatives substituting 



V 



n„ agam 



1 



and the part M without derivatives 



(71) 



4Am ^ I 4At) 



where u = a + 4c + 7, t> = 6 + 4(i — ^ 



The one-loop divergent part of the effective action has the formQ 

= -^tr I (fx ^(]-{Kn-^){Wn){Kn~^){Wn) - {Kn~^){M)] , 

where the matrix product of {Kn~^)^,,^a.(3i (Wn) ^i^^ais, M^^^a/s is supposed. 
A simple calculation gives 



3277% 



J d^x (C4A + C5) 



where 



(72) 



C4 



M 2A + 



( 14A + - + 20 



C5 



5 1 

+ 



/52fc4 V2 8a V 



Since the space-time is fiat the contributions of the Faddeev-Popov ghosts and of the "third" 
ghost are equal to zero. 
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Appendix B 



In this appendix the result of calculation of the two-loop divergent as ^ part of 
the effective action is presented. Also, validity of the general statements of [0, |21| is 
verified. 

Following the algorithm derived in sec. 4b we obtain 



2 (327r%)' 



J d^X (C22i?' + C32i? + C42A + C52) , (73) 



-2 



C22 = ^(200/27a3 - 416/9a2 + 1697/54a - 25/36), 

C32 = ^ {A'(-1/4m2 - 7/2UV - 49/4t;2) 

+A(-3/4M2a"^ - ll/2Mf«"^ - 5Mf + 10/3Ma 
-l/12Ma"^ - 183/20M - - 35v^ + 70/3va 

-IjVlvor^ - 1281/20^;) - Q/lGu^a-^ - l^j^uva-^ 
-3/Suva-^ - l/lQua~^ + 5/2m - bl2v'^a~^ 
-l/lQv'^a-'^ - 2hv'^ + im/?,va - 5/6m~^ - 7/48t;a 
-272/3W - 200/9a2 + 122a - l/24a-2 _ 731/18} , 

C42 = {2A2(u2 + i4^uv + 49t;2) 

+A(6u^a~^ + AAuva'^ + 40ut; + 133/5m + lAv'^a'^ 
+280^;^ + 931/5f ) + 'd/2u^a'^ + 60Mt;a"^ + ^uva''^ 
-Ibua'^ + 5/4ua~^ + 10m + 20^7^0"^ + l/2i;2a"2 
+200i;2 - Sm"^ + 5/12^-2 + 808/3v} , 

C52 = I A(-l/4na-2 _ 5^ _ 7/4^;«-2 _ 35^) 

-15/2ua~^ - 3/8na"^ - b/2va'^ - 5/2i;a~^ - l/8m" 
-50t; + 50/3a - 5/12a-^ + h/Sa'"^ - l/Sa-^ - 79} . 

To show the gauge and parametrization dependence can be absorbed by a field 
renormalization we first remove one-loop divergencies (|56|) by the following redefinition 
of charges and fields 

A-A(1 + 5,A), ^-^(1 + ^1^), 



3 



ai ai{l + 5iai), (3 [3{l + (5i/?), 



where 
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" 327r%' 
24 + C5/2A 



327r% 



327r%ai' 327r%/3' 

and we have introduced a notation Cg for the gauge and parametrization independent 
part of the coefficient C3 

C3 = 4 - J. 

As seen from the above equations renormahzations of the gravitational constant and 
of the metric field can not be separated from each other. This property is inherent to 
any metrical theory of gravity with the Lagrangian containing terms linear in curvature 
and holds at any order of perturbation theory. 

To make the theory finite at the two-loop level we should take into account coun- 
terterms which arise in the second order from the one-loop redefinitions of the charges 
and fields which were made above. Correspondingly, we extract these counterterms 
from the two- loop order result (|7^) rewriting coefficients 032,042 and C52 asQ 

C32 = y 4c4 + ^{A\-1/Au^ - 7/2UV - A9/Av^) + A{-3/Au^a-' 

-n/2uva-^ - 5uv - 133/20M - 7/Av^a'^ - 35v^ - 931/20t;) - 9/16u^a'^ 
-15/2uva-^ - 3/8uva-^ + 15/4m«"^ - S/lGwa^^ - 5/2u - 5/2v^a^^ 
-l/16v^a-^ - 25v^ + 5/Ava-^ - b/ASva'^ - 202/3w - 200/90^ 
+122a - l/2Aa-^ - 731/18}, 

C42 = k^cl/A - 124} + -^{A\u^ + Uuv + 49v^) + A{3u^a-' + 22uva-' 



+20uv + 133/5M + 7v^a'^ + UOv^ + 931/5f) + d/Au'^a''^ + 30uva-^ + 3/2uva 
-Ihua-^ + h/Aua-^ + 10m + l^v'^a'^ + l/Av'^a''^ + lOOt;^ - 5m"^ + 5/12^-^ 
+808/3t; + 400/3^2 - 200a + 5a"^ + l/12a"2 _^ 205/3}, 

C52 = A;2c5(c4 - 44) - -5^{50/3a + h/Aa'^ - Ih/Sa'^ + l/12a-3 + 54}. 
Now it is easy to verify that 



2 



4c32 + C42 - fc^(2c3^C4 + c\/A - 124.) = 
-^{400/9a^ + 288a + 5a"^ - l/12a"^ - 847/9}, 

which means that after subtraction of the counterterms corresponding to the one- 
loop renormalization of charges and fields is made the two-loop divergent part of the 
effective action becomes gauge and parametrization independent on-shell. Therefore 
the gauge and parametrization dependence can be absorbed by a field renormalization 
or by renormalization of the Newtonian constant. 

^^The one-loop redefinitions do not affect the coefficient 022- 
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